Rashba spin orbit interaction in a quantum wire superlattice 
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In this work we study the effects of a longitudinal periodic potential on a parabolic quantum 
wire defined in a two-dimensional electron gas with Rashba spin-orbit interaction. For an infinite 
wire superlattice we find, by direct diagonalization, that the energy gaps are shifted away from 
the usual Bragg planes due to the Rashba spin-orbit interaction. Interestingly, our results show 
that the location of the band gaps in energy can be controlled via the strength of the Rashba 
spin-orbit interaction. We have also calculated the charge conductance through a periodic potential 
of a finite length via the non-equilibrium Green's function method combined with the Landauer 
formalism. We find dips in the conductance that correspond well to the energy gaps of the infinite 
wire superlattice. From the infinite wire energy dispersion, we derive an equation relating the 
location of the conductance dips as a function of the (gate controllable) Fermi energy to the Rashba 
spin-orbit coupling strength. We propose that the strength of the Rashba spin-orbit interaction can 
be extracted via a charge conductance measurement. 

PACS numbers: 71.70.Ej, 85.75.-d, 75.76.+J 



I. INTRODUCTION 

During the last two decades there has been much inter- 
est in using the electron spin in electronic devices. This 
research field, often referred to as spintronics, has already 
made great impact on mctal-bascd information storage 
systems. There are hopes that a similar success can also 
be achieved in semiconductor based systems 1 ^. Manip- 
ulating the spins of the electrons via external magnetic 
fields over nanometer length scales is not considered fea- 
sible. Another, more attractive, method is to use electric 
fields to manipulate electron spins via spin-orbit inter- 
action. The spin-orbit interaction arises from the fact 
that an electron moving in an external electrical field ex- 
periences an effective magnetic field in its own reference 
frame, that in turn couples to its spin via the Zeeman 
effect! 3 -. 

In condensed matter systems, the spin-orbit interac- 
tion is found in crystals with asymmetry in the underly- 
ing structure^. In bulk this is seen in crystals without an 
inversion center (e.g zincblcnde structures) and is termed 
the Dresselhaus spin-orbit interaction 5 -. On the other 
hand the structural asymmetry of the confining poten- 
tial in heterostructures gives rise to the so called Rashba 
term£. The Rashba interaction has practical advantages 
in that it depends on the electronic environment of the 
heterostructure which can be modified in sample fabrica- 
tion and in-situ by gate voltages^. This results in the 
possibility of varying the spin-orbit interaction on the 
nanometer scale. Interestingly, even structurally sym- 
metric heterostructures can present spin-orbit interaction 
provided that coupling between subbands of distinct par- 
ities is allowed^iiS. 

The spin-orbit strength can be measured in a variety of 
different experimental setups^: In a magnetorcsistancc 



measurement via Shubnikov-de Haas oscillationsi&i 2 -"— , 
weak (anti-) localizatio n 13 ^ 5 " —, or electron spin reso- 
nances in semiconducting nanostructures 1 ^— or quan- 
tum dots^ 2 -, or optically via spin relaxation^ 3 -, spin 
precession^, spin-flip Raman scattering 2 ^, or radiation- 
induced magnetoresistance oscillation 
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FIG. 1. Schematic of the proposed experimental setup. The 
top gates produce the parabolic confinement of the quantum 
wire while the finger gates produce the longitudinal periodic 
potential. The Rashba spin-orbit interaction is controlled by 
the backgate. 

In this paper we propose a method for extracting the 
strength of the Rashba spin-orbit interaction via a charge 
conductance measurement. This method does not require 
the use of external magnetic fields or radiation sources. 
We consider a quantum wire modulated by an external 
periodic potential, Fig. [T] Essentially, the method relies 
on the fact that the Rashba-induced shifts of the band 
gap positions in energy dramatically alter the charge con- 
ductance of the superlattice. 

Via direct diagonalization, we determine the band 
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structure of an infinite parabolically confined quantum 
wire. The energy bands clearly show band gaps that 
are renormalized by the Rashba interaction. Interest- 
ingly, the band gaps shift in energy as the strength of the 
Rashba interaction is varied. The location of the band 
gaps are at the crossing of energy bands from adjacent 
Brillouin zones. These energy bands can be calculated 
via an analytical approximation scheme^. 

Using Non-Equilibrium Green's Functions (NEGF) 
and the Landauer formalism we calculate the charge con- 
ductance through a finite region containing both a peri- 
odic potential and the Rashba interaction. In the con- 
ductance we find several dips appearing at different lo- 
cation in energy. Moreover, these conductance dips coin- 
cide with renormalized band gaps of the superlattice. As 
with the band gaps the positions of the conductance dips 
are at the crossing points of the energy bands from next 
neighbor Brillouin zones. For a wide range of Rashba 
spin-orbit interaction strengths some of the conductance 
dips shift linearly in energy as function of the strength of 
the Rashba coupling. For this range we derive a relation, 
Eq. (|I2[) , describing the location of linearly shifting dips 
in energy as a function of Rashba interaction strength. 
The Rashba coupling can therefore be extracted by fit- 
ting the shift of conductance dips via Eq. (fT2"j). Figure 1 
shows a schematic of the proposed experimental setup. 

This paper is organized as follows. In Sec.|TT]we calcu- 
late the energy bands of the infinite wire superlattice via 
direct diagonalization. By using an analytical approxi- 
mation scheme^ we then show how the resulting band 
gaps depend on the Rashba spin-orbit interaction. Via 
the NEGF method and the Landauer formalism we in- 
troduce in Sec. IIIII a numerical scheme to calculate the 
charge conductance through a finite region. This region 
is connected to electron reservoirs to the left and right, 
and contains both a periodic potential and a Rashba spin- 
orbit coupling. In Sec. IIVI we then show that the band 
picture of the infinite wire superlattice, developed in Sec. 
ILT1 is applicable to the finite length periodic potential. We 
close Sec. IIVI with a discussion about possible experimen- 
tal procedures. Lastly, in Sec. [V]we show that our results 
are robust against fluctuations in the strength and width 
of the periodic potential. 



where A is the period of the superlattice. The Hamilto- 
nian that describes this system is 
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were m* is the effective mass, p x and p y are the momen- 
tum operators in the longitudinal and transverse direc- 
tion of the wire, a is the Rashba spin-orbit strength, and 
loq is the confinement frequency of the parabolic poten- 
tial. 

To find the eigenvalues of the Hamiltonian in Eq. ([2]) it 
is convenient to introduce the standard ladder operator 
a of the parabolic confinement and rotate the spin oper- 
ators so that the p x part in the Rashba interaction term 
couples to the a z operator— 
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here a + = a x +ia y is the spin ladder operator, (?r 
is the rescaled Rashba strength and h{k — ^-n) are the 
eigenvalues of the operator p x . In Eq. ([3]) we scale all 
lengths in oscillator length I = ^Jh/m*uiQ and all energies 

in Hloq- We separate the Hamiltonian H into a diagonal 
part, 



Hn 
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and a non-diagonal part 
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II. MODEL SYSTEM: THE INFINITE WIRE 
SUPERLATTICE 

We investigate an infinite quasi-ID parabolic wire with 
a uniform Rashba spin-orbit coupling in the presence of 
a longitudinal modulation described by the potential 

V P (x) = V p0 J2 c ne m ^ x , (1) 

n 



A. Zeroth-order eigenstates and eigenvalues 

The eigenstates of the H Q Hamiltonian are represented 
by the kets \k, m, s). Here m is the quantum number of 
the harmonic transverse energy bands, i.e. the eigenvalue 
of the a' a operator, and s is the eigenvalue of the a z 
operator with s = +1 and s = — 1 denoting the spin up 
and spin down states, respectively. The corresponding 
eigenenergies are 

E^Jk) = i (k - ynj + i + m - sq R (k - ^n^j 
+ c V p0 . (6) 
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FIG. 2. (a) Energy bands of the zeroth-order Hamiltonian 
Ho for I = 50 nm, A = 62 nm= 1.21, ot — 8 meVnm (qjj, = 
0.21/- 1 ), fc* = 2.3r\ c = 1/2, and V p0 = 0.5. (b) Crossing 
of the energy bands corresponding to the |fc, 0, —1), \k, 1, +1), 
and |fc — 2n/X, 0, —1) states, (c) Crossing of the energy bands 
corresponding to the |fc, 1, — 1), \k, 2, +1), \k — 2-n/X, 1, — 1), 
and \k — 2n/X, 0, +1} states, (d) Same as (c) but for higher 
bands, i.e. m — > m + 1. 



A plot of the s (k) energy bands vs k in half of the 
Brillouin zone, i.e., k = . . . n/X, can be seen in Fig. [2] 
(a). 

The Hamiltonian Hi introduces couplings i) between 
the \k,m, s) and the \k, m±l,s) states, due to the 
Rashba interaction, and ii) between the \k,m,s) and 
the \k i n2n / A, m, s) states, due to the periodic poten- 
tial. Note that the coupling is strongest where the energy 
bands corresponding to these states cross each other. For 
some particular Rashba strength a* the energy bands 
associated with the states \k, 0, — 1) and \k, 1, +1) of Hq 
cross at k* = 1/(2^), where = a*m*/h 2 . If we 
choose the period A of the superlattice potential, Eq. ([T]), 
as 



A = 2tt 
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the \k— 27r/A,0, — 1) energy band will also cross at k*. 
This crossing point occurs in energy at 



E * = 7TT + 1 + c oVr 



pO; 



(8) 



see Fig. [2] (b). In the following we refer to this par- 
ticular choice of parameter a* as the reference spin- 
orbit coupling strength. Similar crossing also occur at 
k* with the energy bands associated with the states 
\k,m,—l), \k,m + 1, +1), and \k — 2n/X 7 m, — 1) for en- 
ergies E* + mtvjjQ with the energy band corresponding 
to the state \k — 2ir/X, m — 1, +1) crossing close by. The 
crossings for m — 1 and m — 2 can be seen, respectively, 
in Fig. [1(c) and (d). 



B. Coupling between eigenstates 

We calculate the eigenenergies of the full Hamiltonian, 
H, via direct diagonalization. The resulting energy bands 
are plotted in Fig. [3] (a). The energy bands that cross at 
k* and the one crossing close by, see Sec. Ill A[ corre- 
sponds to the states coupled by H±. A blow up of the 
resulting energy gaps can be seen in Fig.[3](b). At E* the 
coupling results in a double energy gap, see Fig. |3] (c), 
and a triple energy gap at the higher energy crossings, 
see Fig. [3] (d) and (e). In Fig. [3] the spin-orbit strength 
is at the reference value a*. 

Note that for a non-zero Rashba coupling, as in Fig. 
[31 the energy gaps have shifted from the Bragg plane at 
k = tt/X. This is because the spin-orbit interaction shifts 
the wave-number k of the electrons in the longitudinal 
direction and thus renormalizes the locations of interfer- 
ences. 




FIG. 3. (a) Energy bands of the Hamiltonian H = Ho + Hi 
for the same parameters as used in Fig. [5] The energy bands 
are calculated via direct diagonalization of H. (b) A blow up 
focusing on the band gaps of the first three crossings at fee- 
(c) The band gaps formed at the energy band crossing shown 
in Fig.[2](b). In (d) and (e), the energy band crossings shown 
in Fig. [2](c) and (d), respectively. 

When the strength of the Rashba interaction is 
changed the crossing points k c of the energy-bands shift 
and with them the energy gaps. These crossing points 
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can be worked out analytically by using the eigenener- 
gies of the effective Rashba Hamiltonian of the parabolic 
wire2&, which are 



C,t( fc > or) 



(k-nf) 
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and 
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for k - 2tt/A > 0. For k - 2n/\ < e m . s (k,q R ) = 
em,-,(-\k\,qa). In Eq. ® and ([TDJ) 



A m (k,q R ) 
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and we have added the energy constant coVpo result- 
ing from the periodic potential and replaced k with 
k — n2n/X. Note that the energy bands described by Eq. 
© and Eq. (fTO)) are derived for wires without a longitudi- 
nal periodic potential and therefore do not contain energy 
gaps that result from the periodic potential, i.e. they are 
the solution of the H, see Eq. ([3]), with Co 7^ and c„ = 
for n 0. Having determined the crossing points, we in- 
sert them into either of the crossing energy bands, Eq. 
© or Eq. (fTU|) . to obtain the location of the crossing 
points in energy as a function of a. In the appendix we 
present the equation for the crossing point between the 
energy bands i(k,q R ) and sl n _ 1 ^(k,q R ) correspond- 
ing to states |fc,m + 1, +1) and \k — 2ir/\, m, — 1), see 
Eq. (|A3[) . To determine the location in energy of this 
crossing point as a function of a we insert fc c from Eq. 
into either e^^(fc,g R ) or s^^k, q R ). 
Now, as the crossing points shift in energy with a they 
can be thought of as trajectories. The trajectories of the 
crossing points, calculated from Eq. ((9j) and (fT0|) . can 
be seen in Fig. UJ For higher energies there are further 
trajectories. 

We index the trajectories via m and the letters "C", 
"L" , and "R" . The trajectories that we label with the 
letter C only shift a little to the left for a > a* (sec 
lowest horizontal line in Fig. 2]). Relative to the C labeled 
trajectories and again for a > a* , we see trajectories that 
make a large shift to the left and right. Those trajectories 
we, respectively, label with the letters "L" and "R" . 

For spin-orbit strengths close to the reference strength, 
a*, the crossing points follow nonlinear trajectories. As 
the spin-orbit coupling becomes larger or lower than the 
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FIG. 4. (Color online) Trajectories of the crossing points be- 
tween the energy bands. The black dashed lines are the lin- 
earization of the left moving trajectories around ao- We also 
mark the values a* = 8 meVnm and ati = 12 meVnm onto 
the j/-axis. These values are used for the conductance results 
in Fig. El 



reference strength the trajectories quickly become more 
linear. The choice of the reference Rashba spin-orbit 
strength a* determines the period of the periodic po- 
tential, see Eq. Q. 



C. Extracting the Rashba coupling 

We will show in Sec. IIIII that the band gaps appear as 
dips in the charge conductance through a finite periodic 
potential. By fitting the measured energy shift of the con- 
ductance dips to the trajectories in Fig. Slit is possible to 
extract the value of a. The linear parts of the trajectories 
are best suited for fitting. It would therefore be conve- 
nient that the range of a extracted via the fitting is con- 
tained within the linear region. This can be achieved by 
choosing a sufficiently low a*-valuc (and thus A, Eq. [7]). 
We present below a linearized equation for the a- value of 
the crossing point between the energy bands correspond- 
ing to the states \k, m + 1, +1) and \k — 2ir/\, m, —1) as 
a function of Fermi energy. By Taylor expanding around 
some point ao in the linear region of the trajectories and 
making a linear approximation we obtain 



m* I 



mo 



E F - ±V P o - F{m) 
G{m) 



(12) 



here F(m) and G(m) are known functions of the energy 
band index m, see Eq. (|A7|) and Eq. (|A8|) . The derivation 
of Eq. (fT2"]) is shown in the appendix. In Fig.@]we plot as 
dashed lines linearized trajectories described by Eq. (fT2")l 
where we have chosen ao = 20 meVnm. 
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III. FINITE PERIODIC POTENTIAL 

In this section we calculate in the linear response the 
conductance through a finite periodic potential. This is 
done via the Landauer formula together with the NEGF 
method. 



A. The system setup 



a hardwalled wire of width L y with a 



We consider 

transverse parabolic potential. We divide the wire into 
a finite central region of length L x and semi-infinite 
left and right parts, see Fig. [5] The central region in- 
cludes a Rashba spin-orbit interaction described by a 
symmetrized Hamiltonian which is turned on smoothly, 
at both the left and right ends. In the central region we 
also assume a longitudinal periodic potential that repre- 
sents the potential due to the fingergates. We describe 
the central region by the Hamiltonian 

He = ^M+Pl) + \ m *^lv 2 

+ 7^ {<x(x,y),pya x -p x a y } 

1 / /2tt 
+ V p°2 [ 1 - cos \ ~ x 



(13) 



Here {, } denotes an anticommutator. 

The left and right leads contain the same parabolic 
potential as in the center region but neither a Rashba 
spin-orbit interaction nor a periodic potential in the lon- 
gitudinal direction , i.e., they are described by the Hamil- 
tonian 



1 / 2 i 2n,1 * 2 2 



(14) 



The total system, Ht = H\, + He + -Hr, is discretized 
via the finite difference method on a grid of N x x N y 
points with a mesh size a. A schematic of the system can 
be seen in Fig. [5] 



B. The numerical formalism 

Here we use the NEGF method^Zr— to calculate the 
charge conductance. The method requires us to find the 
retarded Green's function of the central region G r c . To 
do this we have to isolate G c from the infinite matrix 
equation describing the retarded Green's function of the 
total system 



(EI - H T )G^ = I. 



(15) 
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FIG. 5. Schematic of the system used in our numerical simu- 
lations. The system is divided into a central region of length 
L x , described by the Hamiltonian He and two semi- infinite 
wires, described by the Hamiltonians H^ and Hr. The total 
system is discretized on a grid with mesh size a. At the left 
and right edges of the central region the Rashba spin-orbit 
coupling is smoothly turned on. 



By separating the total Green's function into a left, right, 
and central part we can write Eq. (|15[) as 



EI - H L -H LC 
— Hql EI — He — Hep 
— H RC EI — H R 

10 0" 
10 
1 



Gf 



G LC G LR 



G 



CR 



Grl G rc Gp 



(16) 



The matrices Hcj = Hjc = tl, couple together the cen- 
tral region to the left (j = L) and right leads (j = R). 
Here t = h/2m*a 2 is the tight-binding hopping parame- 
ter that results from the discretization. Multiplying out 
Eq. (|16p gives nine matrix equations from which we can 
isolate a finite matrix equation for Gjj. This is done 
by treating the contributions from the infinite leads as 
self- energies^— . The matrix equation for G T C is 



(EI-H c -^S j )G r c = I, 



(17) 



and the self-energy of lead j =L, R is Ej 



I f( • G IT 



where Gjj is the Green's function of lead j and is deter- 
mined analytically^. Here He is the discretized Hamil- 
tonian of Eq. ([2) over the central region. All the ma- 
trices are 2N x N y x 2N x N y matrices. In order to save 
computational power only the necessary Green's func- 
tion matrix elements are calculated with the recursive 
Green's function method^. Here we are interested in 
low biases and hence focus on the linear response regime. 
From the Green's function the charge conductance at the 
Fermi Energy Ep can be calculated, via the Fisher-Lee 
relatio: 



,31 



2c 



G(E F ) = — Tr[r R G r c r L G? 



(18) 



where Tj = -2Im(Ej) and G c = (G r c ) T . In what fol- 
lows we use Eq. (fTBj) to calculate the charge conductance 
through our system. 
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C. Numerical values used in the calculation 

In our simulations we consider a Gai-^In^As alloy 
with an effective mass m* = 0.041to o , with m e being 
the bare electron mass. The width of the wire is set as 
L y = 500 nm and its length as L x = 4.92 /im. The oscil- 
lator length is I = 50 nm which corresponds to hoj = 0.74 
meV. For the discretization we use N x x N y = 600 x 63 
points with a = 8.1 nm being the distance between 
nearest neighbors. This corresponds to a tight-binding 
hopping parameter t = 29 meV= 39 Tvjjq. We choose 
the reference Rashba strength as a* = 8 meV nm (i.e. 

= 0.21 which results in an energy band cross- 
ing point at k* = 1/(2^) = 2.3 l^ 1 and a period 
A = tt/ (k* + q R ) = 1.23 I = 7.62 a = 61.5 nm. This cor- 
responds to L x /\ =80 fingergates. The strength of the 
periodic potential is set as V p o = 0.5hxj = 0.37 meV. 

IV. RESULTS 

In Fig. [6] we compare the calculated conductance, Fig. 
[5] (a), through the finite wire to the energy bands of the 
infinite wire, Fig. [5] (b) and (c), as a function of Fermi 
energy. Results for two values of a are presented, a* = 8 
meVnm and ot\ = 12 meVnm. For comparison we also 
plot in Fig. [6] (a) the conductance through a non-periodic 
potential V p (x) = 1/2 V p o for the same values a men- 
tioned above, i.e. the corresponding gapless systems^. 
Each energy gap is associated to its corresponding con- 
ductance dip via a labeled arrow. The labeling on the 
arrows is the same labeling as those of the trajectories 
in Fig. 2J We see that there is a good correspondence 
between the dips in conductance and the gaps in the en- 
ergy band for both a* and ai. This indicates that the 
behavior of the finite periodic potential can be properly 
described with the band model from the infinite periodic 
potential. In Fig. [4] we can sec how the band gaps shift 
in energy as a function of a. We want see if the energy 
shift of conductance dips correlates with the shift of the 
band gaps. This is most easily seen by considering the 
differential conductance 

dG ^ G(E Fi+1 ) - G(E Fi ) 
dE F E-p i+1 — E Fi 

In Fig. [7]we plot as a function of Ep and a. There 
we see the conductance steps of the parabolic wire as rel- 
atively straight vertical trajectories appearing at the tic 
marks on the horizontal axis. The rest of the trajectories 
in Fig. [7] correspond to dips in the charge conductance. 
By comparing the trajectories, that the dips in charge 
conductance form, to the ones of the band gaps in Fig. [4] 
we sec an excellent match. This firmly confirms that the 
conductance dips of the finite wire can be described via 
the band model for the superlattice. We also plot in Fig. 
[7] as light-gray lines the linearized trajectories of crossing 
points between the energy bands corresponding to the 
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FIG. 6. (Color online) Comparison between energy gaps of 
the infinite superlattice and conductance dips of the finite 
periodic potential. In (a) the results from the numerical cal- 
culation of the conductance is shown for ot\ = 12 meVnm, 
and ot* — 8 meVnm. Note that the conductance curve for 
ot\ has been shifted by 2x2e 2 /h for clarity. The results for a 
wire with a non-periodic potential V p (x) = 1/2 V p o is plotted 
with black solid lines. In (b) and (c) the energy bands of the 
infinite wire superlattice are plotted for the same a values as 
in (a), respectively. Note that (c) is the same figure as Fig. 
|3](b), but rotated clockwise by 90°. The correspondence be- 
tween the nine pairs of dips and band gaps are indicated by 
the labeled arrows. 

\k, m + 1, +1) and \k — 27r/A, m, —1) states, see Eq. (fl~2"|l . 
The trajectories are linearized around a=20 meVnm and 
are identical to the straight dashed lines in Fig. @] 

A. Discussion of possible experimental procedures 

From Fig. [7] we can extract the Rashba spin-orbit cou- 
pling by fitting the linear energy shift of the conductance 
dips, via Eq. (|T2"|) . The energy shift, i.e. the change 
in Fermi energy, is controlled by the chemical potential 
of the leads. The backgate that controls the spin-orbit 
strength introduces an extra shift in the Fermi energy 
(due to the electrostatic coupling of the backgate to the 
2DEG). This shift can be compensated by changing the 
chemical potential of the leads by an equal amount. Sim- 
ilar methods have been used to probe the energy spec- 
trum of quantum dots using transport methods^. An- 
other way to compensate for this shift is to use a com- 
bination of back and front gates. This has been experi- 
mentally demonstrated to control the Rashba spin-orbit 
interaction strength without introducing charging in the 
2DEGM. 
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FIG. 7. (Color online) Plot of the differential conductance 
dG/dEp through the finite superlattice as a function of the 
Fermi energy Ep and spin-orbit coupling, a. The linearized 
trajectories of the crosspoints between the energy bands cor- 
responding to the \k, m + 1,+1) and |fc — 2tt/X, m — 1 , +1 } 
states are shown with gray solid lines. The trajectories are 
linearized around ao = 20 meVnm. The Rashba spin-orbit 
values, a* — 8 meVnm and ai = 12 meVnm, correspond- 
ing to the curves in Fig. [6] are marked on the right border. 
Note that 5r = m*a/h 2 l and that the differential conduc- 
tance dG/dEp is in arbitrary units. 



V. VARIATIONS OF PARAMETERS IN THE 
PERIODIC POTENTIAL 

The periodic potential plays a crucial role in in the 
formation and behavior of the energy gaps. We there- 
fore devote this section to study the effect that variations 
of parameters in the periodic potential has on the con- 
ductance. In what follows we consider variations on the 
length and strength of the periodic potential (Sec. IV Aj) 
as well as random fluctuations on the period and strength 
of each finger gate potential (Sec. IV Bp . 

A. Changing the periodic potential parameters 

Here we examine the effects of changing the strength 
of the periodic potential, V p o, and the length L x . In Fig. 
[8] (a) the length of the wire is varied such that it contains 
10, 40, 80, or 160 fingergates with a fixed period A = 63 
nm. There we see the dips easily for > 40 fingergates 
and they become clearer for more fingergates. There is 
not much change between the conductance results for 80 
and 160 fingergates, i.e. w 80 fingergates is an adequate 
number of fingergates. Fig. [8] (b) shows results where 
the strength of the periodic potential has been varied be- 
tween values of 0.2, 0.5, 1.0, and 1.5 hu>o. We see that the 
width of the conductance dips grows with larger V p q. But 



Ep - iVpo [hui ] 

FIG. 8. (Color online) Charge conductance G through the 
finite periodic potential as a function of Fermi energy Ep 
for (a) different values of V^o and (b) different number of 
fingergates. All the conductances were calculated with a — 12 
meVnm. Note that the conductances have been separated by 
2 X 2e 2 /h. 

as Vpo increases the total conductance strength G deteri- 
orates due to interferences. This makes the detection of 
the conductance dips difficult. 

B. Effects of fluctuations in the periodic potential 

The periodic potential used in the previous section can 
be considered as being ideal. A more realistic case would 
be if there were some fluctuations in the potential. To 
verify the robustness of the results in the previous sec- 
tion we rerun the simulations with a non-ideal periodic 
potential. This is done by introducing a 5% Gaussian 
error in both the height and length of each potential hill 
created by the fingergates. Figure [§](b) shows that extra 
noise is added to our conductance result; still the tra- 
jectories of the conductance dips are fairly visible. To 
counter against the noise introduced by fluctuations we 
could add extra periods to the periodic potential. This 
helps averaging out the noise, as can be seen in Fig. [9] (a) 
where we have quadrupled the length of the wire and the 
number of fingergates. This (self) averaging is however 
slow and requires a large number of fingergates. Another 
solution, as we are mainly interested in the slope of the 
crosspoint trajectories, is actually to reduce the number 
of fingergates to an optimal number. As can be seen 
in Fig. [5] (c) we obtain many extra trajectories resulting 
from the imprecision of the period. But the slope of the 
trajectories of the conductance dips can be easily fitted 
for a > 10. As expected if we go down to as few as 10 
fingergates the effects of the periodic potential is nearly 
washed out. 
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FIG. 9. Plots of the differential conductance dG/dEp through 
the finite periodic potential as a function of Fermi energy 
Ep and spin-orbit coupling, a. Four cases are considered: A 
system with (a) 320 , (b) 80, (c) 40, and (d) 10 fingergates. In 
all the plots we introduce a 5% Gaussian error in the length 
and height of the potential in each period. 



VI. CONCLUSION 

We have studied a parabolic quantum wire with 
Rashba spin-orbit interaction and a longitudinal periodic 
potential. Wc find that the energy gaps resulting from 
the periodic potential split up and shift from the Bragg 
plane due to the Rashba spin-orbit interaction. Above a 
certain spin-orbit strength some of the new energy gaps 
shift linearly in energy as a function of the Rashba spin- 
orbit strength. We propose that this effect can be used to 
measure the change in the Rashba spin-orbit strength, e. 
g. resulting from a voltage gate^. The energy gaps result 
in dips in the charge conductance. The energy shift of 
these dips can be fitted, by an analytical equation that 
we derive, Eq. (|T2"j) . to extract the strength of the Rashba 
spin-orbit interaction. The advantage of the this method 
is that it only requires conductance measurement and not 
any external magnetic or radiation source. 



Appendix A: Energy crossing points 

In this appendix we calculate the crossing point k c of 
the |/c,to, — 1) and \k — 2tt/\, m — 1, +1) energy bands. 



n -1 2tt 27T 2 

£ m ,li k > <?R) - e m-l,f( fc ) OR) = ~ k - — + 1 



A 2 



-A(k, q R ) - A m+1 (fc, q R ) + A m {-k + 2-, q R ) = 



(Al) 



for k. Here we used Eq. © and (|T0|) and have defined 

1 1 



A(k,q R ) = f 

2 V I- 2q R k 



— QR 



1 + 2q R k 



(A2) 



Equation IA1I is a nonlinear equation which is hard to 
solve directly. An easier way is to linearize the A(k) 
and A m (k) functions around the point k = k*. This 
approximation introduces little errors and allows us to 
write the crossing point as 



_ - 1 - B(m, g R ) + C(m, g R )k* 



where 



B(m,q R ) 



f + C(m,q R ) 



■ A{k*,q R ) - A m+1 (k*,q R ) 
A m (-k* + 2n/\,q R ), 



(A3) 



(A4) 



and 



C(m,q R ) 



dA(k,q R ) 



dA m+1 (k,q R ) 



Ok 

dA m (k,q R ) 



Ok 



Ok 



(A5) 



k— — fcc+2-7r/A 



The trajectories of the crossing points across the 
energy-Rashba spin-orbit coupling surface are then 
£ m,i( fc c, <7r) or £^ t _ 1)t (fc c ,?R). The crossing points and 
their trajectories for the other energy bands, namely 
the ones corresponding to the states \k,m + 1, +1) and 
\k — 2ir/\, m, —1), the \k, m, —1) and \k — 2n/X, m, — 1), 
and the \k, m + 1, +1) and \k— 2-K/\ 1 m— 1, +1), can be 
worked out in the same way. 

The k value of the crossing point between the 
£ m j,(^> Qr) an d e m-i 1r) t> an ds remains almost con- 
stant as a function of a. This is because these bands 
move at similar rate in opposite directions in k space as 
a is changed. We take an advantage of this when finding 
a linear equation around some point «o in the straight 
segments of the trajectories. We know that k* is a cross- 
ing point for a = a* so we can, instead of using Eq. 
(|A3|) . make the approximation that k c w k* for all a. We 
then linearize ^(k* ,q R ) by Taylor expanding around 
ctQ. Solving for q R in ^(k*,q R ) and scaling back to a 
we obtain 



TO* I 



rO 



F(m) 



G(m) 



(A6) 
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where 



F(m) 



\iK) 



1 - 



<2ro 



2(1 + 2q R0 k* 



3.7 + m 

^6.0 x 10- 5 m 3 - 4.1 x 10- 



- A 



m+l 



0.20m- 



-0.79, 
(A7) 



In the last steps of the Eq. (|A7[) and Eq. (|A8|) we have 
plugged in the values k* = 2.35 and (?ro = 0.537, which 
correspond to a* = 8 meVnm and ao — 20 meVnm. 
Note that all lengths are scaled in the oscillator length 
/ = ^/h/m*LUo and all energies in hu>a. 
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